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Abstract: We study methods for reconstructing the momenta of invisible particles 
in cascade decay chains at hadron colliders. We focus on scenarios, such as SUSY and 
UED, in which new physics particles are pair produced. Their subsequent decays lead 
to two decay chains ending with neutral stable particles escaping detection. Assuming 
that the masses of the decaying particles are already measured, we obtain the momenta 
by imposing the mass-shell constraints. Using this information, we develop techniques 
of determining spins of particles in theories beyond the standard model. Unlike the 
methods relying on Lorentz invariant variables, this method can be used to determine 
the spin of the particle which initiates the decay chain. We present two complementary 
ways of applying our method by using more inclusive variables relying on kinematic 
information from one decay chain, as well as constructing correlation variables based 
on the kinematics of both decay chains in the same event. 
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1. Introduction 

The operation of the Large Hadron Collider (LHC) starts a new era in high energy 
physics of direct exploration into the TeV scale. New physics beyond the Standard 
Model (SM) is strongly expected to occur at the TeV scale because of the hierarchy 
problem. Another strong hint for new physics at the TeV scale comes from the dark 
matter (DM) in the universe. It is now well-established that ~ 23% of the total energy 
of the whole universe is made of dark matter, and it cannot be accounted for from any 
SM particles. The leading candidate for the dark matter is a new weakly interacting 
massive particle (WIMP) with a mass in the range of ~ 10 GeV to a few TeV^. The 
thermal relic of such a particle from the Big Bang can give the right amount of dark 
matter in the universe if its interactions with SM particles and itself are of the similar 

^See Ref. Q for a recent review. 
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strength of the weak interaction. To be stable it should be charged under a new 
symmetry {e.g., a Z2 parity as the simplest example). 

Many scenarios of TeV new physics beyond the Standard Model have been proposed 
to address the hierarchy problem. They often also contain a dark matter candidate, 
e.g., supersymmetry (SUSY) with i?-parity, Universal Extra Dimensions (UEDs) with 
Kaluza Klein (KK) parity [Q, ^, ^, little Higgs models with T-parity warped extra 
dimensions with a symmetry and so on. A common feature of these models is that 
there are other new particles charged under the same new symmetry which protects 
the stability of the DM particle. These new particles may be pair produced copiously 
at the LHC. After production they will go through cascade decays to the lightest one 
which escapes the detector. Therefore, many of these different models can give rise to 
similar collider signatures, i.e., jets/leptons with missing transverse momentum. It is 
important to be able to distinguish different models if such experimental signals are 
found. In particular, the spin measurements are essential to distinguish SUSY where 
the spins of the SM particles and their superpartners differ by 1/2, from other models 
where the partners have the same spins as the corresponding SM particles. 

To determine the spin of a particle in a more model-independent way, we need to 
examine the angular distributions of its production or decay. Although theoretically 
well-motivated, these collider signatures with missing transverse momentum pose a 
serious challenge to such experimental measurements at hadron colliders. Because the 
new particles are pair-produced, there are at least two missing particles (one from 
each chain) in every event. For any given single event there is not enough information 
to reconstruct the full kinematics without additional information. In a long decay 
chain, the polar angle of the decay of the intermediate particle in its rest frame is 
directly related to the Lorentz invariant mass combination of the visible particles of the 
decay chain. One can use the invariant mass distribution to determine the spin of the 
intermediate particle without fully reconstructing the kinematics if certain conditions 
are satisfied. Many of the spin determination methods in the literature are based 
on this observation @, |, [n|, |13|, |15|, |16l 0, |T|, ig, |2|. However, it 



cannot be used to completely determine the spins of the first and the last particles in 
a decay chain directly. If the momenta of the invisible particles of each event can be 
reconstructed, then one can boost the event to any frame and examine any relevant 
kinematic distributions. The spin of a particle can be determined from the azimuthal- 
angle correlations as well as the polar angle distributions. In particular, one can 
determine the spin of the first particle in a decay chain by looking at the angular 
distributions of its production or decay. 

It is possible to reconstruct the momenta of the invisible particles if there are 
enough constraints to match the number of unknown kinematic variables in the event. 
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e.g., if there are enough mass shell constraints and the masses of the particles in the 
decay chains are already known. Measuring the masses of the particles in a decay 
chain with missing transverse momentum itself is a non-trivial task, as there is no 
invariant mass peak and the visible momenta are more sensitive to the mass differences 
than the absolute masses. There have been many research efforts recently in mass 
determinations for various event topologies. A lot of progress has been made and many 
new methods have been proposed based on various kinematic variables and constraints 



(for a review, see ||23|). We expect that the masses of the new particles can be quite 
accurately determined if a substantial clean signal sample can be isolated and the visible 
momenta are well measured. In particular, for extended decay chains we consider in 
this article, the masses can be determined with a few percent errors using a few hundred 
events [|2|, ^, ^. 

The rest of this paper is organized as follows. In Sec. |^, we discuss single-chain 
and double-chain event topologies which can be kinematically reconstructed with the 
mass measurements. In Sec. ^ we discuss spin determination from a single decay 
chain. In particular, we identify a case where we can determine the spin of the first 
decaying particle, which cannot be measured by the invariant mass technique. In 
Sec. we discuss the double-chain techniques, where additional information from the 
spin correlation between the two decay chains may be used. Further discussion and 
conclusions are drawn in Sec. |[ 



2. General Considerations of Event Reconstruction 

As we discussed in the Introduction, it is very useful if we can reconstruct the invisible 
particles' momenta from the visible momenta and available kinematic constraints. In 
this section, we give a general counting of constraints for different topologies of events 
with missing transverse momentum, and discuss the corresponding methods for event 
reconstruction. As we will see, this depends on whether we are examining a single decay 
chain or both decay chains in an event, and whether the system is under-constrained, 
exactly-solvable or over-constrained. 

The unknowns in the problem are the 4-momenta of the missing particles. Assum- 
ing that there is only one missing particle in each decay chain, we have 4 unknowns 
for each event if we want to reconstruct only one of the decay chains, and 8 unknowns 
if we want to reconstruct both decay chains. As we mentioned in the Introduction, 
we will assume that the masses of all new particles in the decay chains are already 
measured with some errors. Each on-shell particle then contributes a constraint on the 
missing momenta. These "mass-shell constraints" are available for both single-chain 
and double-chain cases. For the double-chain case, two more constraints are available 
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from the measured missing transverse momentum if there are no extra missing parti- 
cles. It is then straightforward to count the number of constraints needed for event 
reconstruction. For the single-chain case, we need 4 mass shell constraints to solve the 
system. (We say the system is exactly-solvable when the number of unknowns is equal 
to the number of constraints.) This corresponds to a decay chain with 3 visible particles 
(including particles decaying further but not introducing extra invisible particles, such 
as a Z-boson decaying to charged leptons/quarks) if all decays are two-body. If the 
decay chain is longer, we have an over-constrained system and we can employ a likeli- 
hood method to obtain the best-fit missing momenta. If the decay chain is shorter, the 
missing momenta cannot be fully reconstructed. Similarly, in the double-chain case, we 
need 6 mass-shell constraints, which, together with the constraints from the measured 
missing transverse momentum, allow us to solve the system. An example is the case 
with two on-shell decays for each decay chain. This occurs for ti pair production in the 
dilepton decay channel and can be used for determining the spin of a ti resonance ||28|| . 
Again, if the decay chains are longer (shorter), we have an over- (under-) constrained 
system. 

The single-chain case and double-chain case also differ in the available spin cor- 
relation information. For the single-chain case, the relevant quantity is the angular 
distribution of the decay products. Here, the decaying particle can be the first particle 
or any of the intermediate particles in the decay chain. In order to have a non-uniform 
angular distribution for the decay products, the decaying particle needs be polarized. 
In addition, if the decaying particle has spin 1/2, the coupling responsible for the decay 
needs be chiral. Of course, this information is also available in the double-chain case. 
On the other hand, if both decay chains are reconstructed, we obtain extra informa- 
tion unavailable in the single-chain case, namely, the spin correlations between the two 
decay chains. 

In this paper, we focus on the cases where the system is exactly-solvable or over- 
constrained^. In particular, we analyze in detail decay chains with three visible par- 
ticles, in both the single-chain case and the double-chain case. We will discuss the 
corresponding event reconstruction methods and related issues. It is straightforward 
to generalize the methods to other event topologies. 

3. Single Chain Technique 

3.1 Angular distribution of decay products 

A decay chain with 3 on-shell decays is shown in Fig. |1]. The particles Z (not to be 

^For an under-constrained system, although the missing momenta cannot be fully reconstructed, 
correlations between spin and kinematics often exist, see, for example, Ref. |3T[| . 
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Figure 1: A decay chain with 3 visible SM particles. The final-state particles are labeled 1 
through 4 with 1 denoting the missing particle and 2, 3, 4 denoting visible SM particles. The 
new particles are called Z, Y, X and and assumed on shell. 



confused with the Z boson, which we denote Z^), Y, X and are assumed to be 
on-shell with masses mz, my, rrix and rriN. Before describing the details of event 
reconstruction for this topology, we first discuss how to observe the spin correlation in 
the single-chain case once the missing particle's momentum is obtained, and compare 
it with the invariant mass method studied in the literature. 

The basic idea for observing spin correlation from a particle decay is as follows: 
suppose the decaying particle is polarized and the coupling responsible for the decay is 
chiral, then the decay products will have a non-uniform angular distribution in the rest 
frame of the decaying particle. The two daughter particles' momenta are back to back 
in the rest frame of the mother particle and we use 6 = 6{m, d) to denote the angle 
between either of them and the polarization axis of the mother particle. Here m denotes 
the mother particle and d denotes the daughter particle. The probability density of the 
decay is a polynomial in cos 9 of order 25", where S is the spin of the mother particle. 
Note that the coefficients of the polynomial depend on the spin density matrix of the 
decaying particle, the coupling responsible for the decay and also the axis one chooses 
to evaluate the angle 9. In special cases, for example, when a fermion decays through 
a vector-like coupling, the coefficient (s) of the leading order term(s) could be vanishing 
or too small, giving a polynomial of order lower than 25*. Therefore, using this method 
we can only set a lower bound on the decaying particle's spin. When the polarization 
axis is coincident with the direction of the mother particle's initial momentum, we say 
that the particle is polarized in the helicity basis, and denote the angle defined above 
by 9hei{m, d). With full reconstruction of the kinematics of the event, 9hei{fn, d) can be 
simply obtained by boosting to the rest frame of the mother particle. 

One may get a polarized particle if itself comes from the decay of another particle 
through a chiral vertex. For example, we consider the particle Y in Fig. |l|, which comes 
from a two-body decay of the particle Z. It then decays to two particles, X and 3. 
In the rest frame of Z, particle Y and particle 4 move in opposite directions and Y is 
polarized along that direction. Now we can boost the system to the rest frame of Y 
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Figure 2: The decay Z — )■ y — t- X viewed in the Y rest frame. 



(Fig. 1^). In this frame one can see that the angle between particles 3 and 4 is simply 
6*34 = — 9hei{y, 3) since the direction of the particle 4 is unchanged under the boost. 
The combined invariant mass of visible particles 4 and 3 can be easily calculated in 
this frame and it is related to the angle 9hei{y, 3) by 

ml, = {Pz+PA? = limT^'fa - cos ^34) = limT^'fa + coseUY.^)). (3.1) 



where 



1 1 6y 



This fact has been used for spin measurements in Ref. [|, |, |1^, 0, |1^, 0, [T^, [16 



1^, |T8|, [T9|, [2^, ^ . The advantage of the invariant mass method is that the distribution 
can be obtained without complete event reconstruction. However, it is also clear that it 
requires the particle to come from a heavier particle decay, hence it can only be applied 
to the spin determinations of the intermediate particles directly in a decay chain^. While 
the last particle never decays and it is hard to obtain its spin information directly, the 
spin of the first particle may be determined directly if its polarization already exists at 
the production level. There is no experimentally measurable Lorentz-invariant quantity 
related to the angle in which we are interested in this case. On the other hand, by 
reconstructing the missing particle's momentum, we can directly examine the angle 
Ohei{m,d) of the first particle decay in a decay chain and extract its spin information, 
which we describe in the following subsections. 

3.2 Momentum reconstruction 

We now describe the momentum reconstruction for a decay chain as shown in Fig. |I|. 



■^Of course once the spins of the intermediate particles are known, we can determine whether the 
first (last) particle is a boson or fermion. 



- 6 - 



The mass shell constraints give the following equations: 



{Pl+P2f 
{P1+P2+P3Y 




(3.3) 



where pi is the four-momentum of the invisible particle, P2,P3,P4 are the four-momenta 
of the visible SM particles, m^, mx, my, mz are the measured masses. It is easy to see 
by taking the differences that these equations can be simplified to 3 linear equations 
plus a quadratic equation for the invisible momentum. Therefore, the system of these 
equations always admits two solutions, with the number of real solutions being or 
2. The solutions can become complex if we use wrong combinations of the visible 
particles or the experimental smearing is too large. Such "bad" events or combinations 
should be eliminated by requiring the solutions to be real. However, in the classes 
of models under consideration, there is always another decay chain beside the one we 
wish to study, which can also contain similar final state particles. Assigning a final 
state particle to the "wrong" decay chain sometimes can also yield real solutions. In 
practice, separating out such contaminations can be very challenging. In our study, 
we choose to accept all real solutions and add them with equal weight. As shown in 
our case study below, making such a choice does not prevent us from extracting spin 
information. A more careful treatment of such combinatorial contamination should be 
able to further enhance the spin differentiation power. 

3.3 Application to decay chains with photons 

For an illustration, we apply the above method on a decay chain in a gauge medi- 
ated SUSY breaking (GMSB) model and its UED counterpart. These models are 
chosen simply because they give rise to the event topologies which are suitable for 
the spin studies with our method. GMSB is characterized by a gravitino LSP (light- 
est supersymmetric particle). We assume that the gravitino has a mass ~eV which 
is essentially massless compared with the detector resolution. The next lightest su- 
persymmetric particle (NLSP) is assumed to be a Bino-like neutralino which decays 
promptly to a photon and the gravitino. We are interested in the following decay chain: 
xl iniR ^r^rX^ ^R^RlG (Fig. | (a)), where X2, Xi and G are respectively 
the second neutralino, the right-handed slepton, the lightest neutralino and the grav- 
itino, corresponding to the particles Z, Y , X and in Fig. |l]. We set the mass of the 
gravitino to be zero and choose the other particles' masses and interactions according 
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Figure 3: Decay chains with photons in GMSB and UED. 



to the Minimal Supersymmetric Standard Model (MSSM) point SPSla"^ Thus the 
masses of the particles Z, F, X and N are 181, 143 , 97 and GeV. As denoted by 
the subscripts, X2 o^ly decays to right-handed leptons. Therefore, if X2 is polarized, 
we will be able to determine its spin. 

For comparison, we consider a similar decay chain in UED, Z^^"* — )■ ^l^l"* — ^ 
^l^lT*-^-* — ^ ^L^Ll4', where , 7^^^ denote the first KK modes of the gauge 

bosons and leptons, and is a scalar field and the lightest KK-odd particle (LKP). 
Note that is absent in the minimal 5-dimensional (5D) UED model, but exists in 
6D UED models as a scalar KK partner of the hypercharge gauge boson (dubbed 
B^j^'^^ in Ref. ||35|). A scalar LKP also exists in an extension of the 5D UED model with 
an additional gauged Peccei-Quinn (PQ) U{1) symmetry as the zero mode scalar 



partner of the PQ gauge boson (denoted B^). If the NLKP (next lightest KK-odd 
particle) is the KK-photon, it decays to a photon (or Z-boson) and B^. Since that B^ 
can be very light (< GeV) while B^j^'^^ in the 6D UED model has to be massive (> 
lOOGeV), the PQ-UED model can mimic more closely the signature of GMSB than the 
6D UED model. For our purpose, the model subtleties are unimportant and we only 
need to fix the spins according to the model and specify the masses and couplings of the 



"'SPSla is actually a model point of the minimal supergravity (mSUGRA) instead of GMSB. Note 
also that Tevatron results have put stringent constraints on events with photons plus missing transverse 
momentum js^ , essentially ruling out a spectrum with such low masses. We just use this spectrum 
for illustration because of its clean chiral structure. It is straightforward to apply the same method 
on other model points. 
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particles. We will choose the masses to be the same as the SUSY case. For simplicity, 
we assume zji^ is purely Wji^^^ so that its couplings to fermions are purely left-handed: 



CDqY^^q^'^ZW+h.c. (3.4) 
The KK-photon decays to (p through a vertex in the form 

where c is a coupling constant. 

The process under consideration is the neutralino/chargino (x^/xf) pair produc- 
tion in GMSB and the UED counterpart KK-Z^/KK-VT^ (zj:^^ /W^t^^^) pair produc- 
tion. Thus X2/ is the first particle in the decay chain. In the lab frame, X2 more 
left-handed. This can be seen as follows: we consider the process ud — )■ x^Xi", the 
other process ud — )■ X2X1 is similar. The process is dominated by the s-channel W^^- 
exchange diagram so that the initial u{d) is left-handed (right-handed). Therefore, in 
the center-of-mass frame, X2 more left-handed (right-handed) in the forward (back- 
ward) direction with respect to the u quark due to angular momentum conservation. 
In the center-of-mass frame, X2 equally left-handed and right-handed^. However, the 
system tends to have a large boost along the u quark direction, changing some of the 
right-handed X2 to left-handed. Therefore we have more left-handed X2 than right- 
handed in the lab frame. Similarly, in the UED case, KK-Z^ is negatively polarized in 
the lab frame. 

Wrong combinations involving particles in the other decay chain can also yield 
real solutions and contaminate the distributions. For SPSla, 95% of xf decays to 
a stau (r) and a neutrino. Therefore, we assume the chargino decays according to 
xf — — )■ Tv^Xi ^ Tv^^G. Accordingly, we let KK-PF^ decay through a KK-r to 
the KK-photon, which then decays to the scalar (p and a photon. Therefore, each event 
contains two photons, two opposite-sign same-flavor leptons and a (hadronic) r, which 
amount to a 4-fold ambiguity for assigning the positions of the two photons and the two 
leptons. We will give a more detailed assessment of this combinatorial contamination 



in Sec. 3.4. 



The events are generated with Herwig++ 2.4.2 at the parton level^ for 14 TeV 
pp collision. For simplicity, we have turned off initial/final state radiations. The final 

^This is not exactly true — X2 is slightly polarized in the center-of-mass frame due to its small 
Higgsino component, see Ref. 

^Only the minimal UED is available in the official Herwig++ code. We have adjusted the code to 
allow a generic particle spectrum. We have also added the scalar field 4> and let the KK-photon decay 



through the couphng Eq. (B.5). The decay Xi — > jG is performed according to phase space since the 



decay products have a uniform angular distribution. 



- 9 - 



state radiation is small since the final state particles are either leptons or photons. The 
effect of the initial state radiation is to give the whole ^/Xi or ^^-Z^/KK-W^ system 
a boost, which do not change qualitatively any of the results presented below. 
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Figure 4: Angle between the near lepton and ^/ Z^^ from Monte Carlo truth. Left: cos 9^+ 
for UED and cosO^- for SUSY; right: cosO^- for UED and cos0£+ for SUSY. The number of 
events is normalized to 10k in 20 bins for all histograms. 



We are interested in the spin of the first particle in the decay chain, i.e., x^/ Z\i . 
Therefore, we examine the angle On = 9hei{Z, i) as discussed above, where i is the "near" 
lepton, namely, the lepton directly from the X2/'^m^'' two-body decay. In Fig. |, we 
show the distributions for cos 6c from the Monte Carlo. For completeness, we draw the 
distributions separately for positive and negative near leptons . Because the coupling 
in the SUSY case is right-handed while in the UED case is left-handed, the distribution 
of 9i+ {6i-) in SUSY should be compared to that of 6^- {Oi+) in UED, which we put in 
the same figure. It is clear that the distribution is linear for SUSY and quadratic for 
UED, corresponding to spin-1/2 and spin-1 particles. We have normalized the number 
of events for each distribution to 10k (20 bins), although we have used more events to 
produce the smooth distributions. Note that it is unnecessary to make this distinction 
based on lepton charge, if we limit our goal to differentiate SUSY from UED, since we 
only need to distinguish linear vs quadratic behavior in this case. On the other hand, 
the slope does carry the information of chirality of the coupling in this case. 

We first apply the event reconstruction method on events without any experimental 
cuts or smearing. The distributions including all solutions with equal weight are shown 
in Fig. ^ (a). For comparison, we have also performed the reconstruction for events 
with the same mass spectrum and 2 — )■ 2 differential cross-section (for UED and SUSY 
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Figure 5: Reconstructed cosOi. Exact momenta without experimental smearing are used. 
All combinations included. Only showing cos0^+ for UED and cosO^- for SUSY. The UED 
(SUSY) PS distribution is obtained using the UED (SUSY) 2 — )• 2 differential cross-section 
with all decays performed according to phase space. Left: before subtracting PS distribu- 
tions; right: after subtracting PS distributions. We normalize the number of solutions for all 
histograms on the left panel to 10k, and do the subtractions to obtain the distributions on 
the right panel without further normalization. 



respectively), but with all particles decayed according to phase space (PS). Comparing 
Fig. |^(a) and Fig. ^, one can see that the distributions have been distorted from the 
Monte Carlo truth by wrong combinations and wrong solutions, and the theoretical 
linear and quadratic functions of Fig. ^ are lost. Nevertheless, the distributions of UED, 
SUSY and PS are clearly distinguishable. We can also retrieve (some of) the theoretical 
behavior by subtracting the UED and SUSY distributions from the corresponding PS 
ones, which are shown in Fig. ^ (b). The subtracted distributions are much closer to 
the original ones, although the contamination cannot be completely removed. Note 
that SUSY and UED give rise to different PS distributions, which can be attributed to 
the difference in the differential production cross section, as discussed in more detail 
later in Sec. p.4| . A further potential obstacle is that, in practice, we do not know which 
PS distribution to compare to. However, as shown in Fig. ^(b), the distributions are 
still distinguishable even if we made the wrong subtraction. 

The actual distribution observed in a collider detector is also subject to modifi- 
cations from experimental smearing, cuts, efficiency, etc. We simulate the detector 
response using a simplified approximation described in the Appendix, taking into ac- 
count the detector coverage and momentum resolution^. The cuts on px (> 10 GeV 



^Of course, our detector simulation is far from a complete one, which has to include effects such as 
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Figure 6: Reconstructed cosOi. Experimental smearing applied. All combinations included. 
Only showing cos0^+ for UED and cos0^- for SUSY. Left: before subtracting PS distribu- 
tion; right: after subtracting PS distribution. We normalize the number of solutions for all 
histograms on the left panel to 10k, and do the subtractions to obtain the distributions on 
the right panel without further normalization. 



for both leptons and photons) and \ri\ (< 2.4 for leptons and < 3.0 for photons) reduce 
the number of events to 82% for UED and 63% for SUSY. The UED efficiency is larger 
because the particles have higher pt, as explained later in Sec. P^ . Since the visible 
particles in our example are either leptons or photons, both of which have good reso- 
lutions in a collider detector, the experimentally smeared distributions (Fig. P) are not 
significantly different from those using exact momenta. 

When producing Fig. H, we have used the correct masses to obtain the solutions. 
In practice, the masses are measured with errors which could alter the distributions if 
they are significantly different from the correct values. However, for decay chains with 
multiple leptons and photons as we are considering, we expect good resolutions of mass 
measurements. To estimate the effect of mass measurement errors, we shift the input 
masses by +5 GeV and repeat the above procedure. The resulting distributions are 
given in Fig. ^, showing only tiny shifts from the distributions in Fig. More impor- 
tantly, it shows that such errors in mass measurements do not change the distinction 
between the SUSY and UED distributions. 

3.4 Combinatorics and more global information 

It is interesting to examine the number of real solutions decomposed according to dif- 
ferent combinations, which is shown in Table ||. As mentioned before, we can eliminate 

trigger efficiency, niis-identification rate, isolation cuts, etc. 
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(a) (b) 

Figure 7: The same plots as in Fig. 6, but with ah input masses shifted by +5 GeV from the 
correct values. The original distributions in Fig. 6 (b) are superimposed in the right panel 
for comparison. 





correct 


wrong photon 


wrong lepton 


wrong photon & wrong lepton 


UED 


2 


0.13 


1.21 


0.12 


SUSY 


2 


0.70 


1.20 


0.67 


UED PS 


2 


0.13 


1.21 


0.12 


SUSY PS 


2 


0.70 


1.21 


0.67 



Table 1: Average number of solutions per event for various combinations. Exact momenta 
without smearing are used. The correct combination always yields two solutions. "Wrong 
photon" means that we have used the photon from the wrong decay chain. "Wrong lepton" 
means that the two leptons in the decay chain are interchanged from their correct positions. 



some wrong combinations by requiring the solutions to be real. Nevertheless, there are 
still significant contributions. From Table |l|, we see that the number of solutions from 
wrong lepton combinations is similar for all cases. On the other hand, the number of 
solutions from wrong photons is sensitive to the kinematics of the other decay chain, in 
this case, the pt of the wrong photon (Fig. H). From the right panel of Fig. ||, we see 
that the pr distribution of the wrong photon is almost identical for SUSY and SUSY 
PS (the same is true for UED and UED PS), but very different between SUSY and 
UED. This is due to the distinction in the 2 — )■ 2 differential cross-section for SUSY 
and UED. As shown in the left panel of Fig. ^ the productions in UED events tend to 
be more central, leading to a harder photon spectrum. 

It may be possible to develop more sophisticated methods to reduce the wrong 
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300 



Figure 8: Left: Production angle in the neutralino chargino center of mass frame. Right: 
PT distributions of the wrong photons. The UED PS (SUSY PS) distributions are shifted 
horizontally by 5 GeV to distinguish from the UED (SUSY) distribution. 



combinations using additional information of the other decay chain. To achieve that 
more precise knowledge of all possible decay chains, such as masses and couplings of all 
particles involved, are often necessary. However, the strategy would be highly model- 
dependent in order to be more effective. Although it could be a useful step in practice, 
we will not pursue this further complication in our analysis. 

Finally, we note that the difference in production angle shown in the left panel of 
Fig. ^ is directly correlated with the difference in spin of the particles. Hence, this 
production angle itself is a very good variable for spin measurement. However, in this 
particular channel, this angle is not readily reconstructable due to the presence of a 
neutrino in the chargino decay chain. Although less direct, as shown above, the pt 
distribution of the wrong photon is sensitive to the production angle, and it can be 
used for spin measurement. In the next section, we will explore cases where we can 
reconstruct both decay chains, then the production angle provides a direct probe of the 
spin. 



4. Double Chain Techniques 

In this section, we present techniques based on reconstruction of both of the decay 
chains. In principle, they are applicable to any event topology with enough constraints 
to solve the kinematics of both decay chains. As a demonstration, we focus on the case 
with two identical decay chains, as shown in Fig. Event reconstruction has been 
considered for this event topology in Refs. p6| , The purpose there is to determine 
the masses of the particles in the decay chain, while the goal here is to find the best-fit 
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8 6 4 

Figure 9: An event with two decay chains, each containing 3 visible SM particles. The 
final-state particles are labeled 1 through 8 with 1 and 2 denoting the two missing particles 
and 3-8 denoting visible SM particles. The new particles are called Z, Y, X and A'^ and 
assumed on-shell. 



momenta for the invisible particles assuming all masses are known (with uncertainties). 
The invisible particles' momenta are obtained as follows. First, we have 8 equations 
from the mass-shell constraints of the 8 on-shell particles in the two decay chains. In 
addition, if the only missing particles in the events are the two neutral particles at the 
end of the decay chains, we have two additional constraints, 

Pi+P2=f, V\+Vl=f- (4.1) 

Therefore, we have 10 equations and 8 unknowns and the system is over-constrained. 
Using the uncertainties in mass measurement given in Refs. [^, together with the 



experimental errors for the visible momenta, we perform a likelihood fit to find event 
by event the best fit momenta of the missing particles. We describe the fitting procedure 
in Appendix 0. 

Once the momenta of the missing particles are reconstructed, we can of course 
obtain the angular distributions of the decay products as in the single chain case. More 
interestingly, we also obtain information unavailable in the single chain case, which we 
illustrate by applying the method on sbottom/KK-bottom pair productions. 

4.1 Application to sbottom/KK-bottom pair productions 

In this process, each event contains two sbottoms/KK-bottoms. The sbottom is as- 
sumed to decay in the following decay chain: 

h^hxl-^Ul^hU^,, (4.2) 
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and a similar decay chain occurs for the KK-bottom. Notice that we have enough 
constraints here to carry out a single chain analysis, using the method presented in 
Sec. 1^. However, in this case, the single chain analysis will not reveal the spin of the 
sbottom or KK-bottom. The decay product of sbottom will give a flat distribution in 
Ohei- At the same time, the KK-bottoms are produced mostly through their couplings 
to gluon which is vector-like. Therefore, the KK-bottom is almost all unpolarized and 
the 9hei distribution of the decay products is fiat as well. Therefore, in this case, we can 
only get the spin information of the sbottom/KK-bottom from a double chain analysis. 

We use the SUSY particle spectrum of SPSla, and we set the UED mass spectrum 
to be the same. The masses for particles Y, X and A^, corresponding to 6, X21 ^1 ^^^1 
X? (or similar KK states), are then {515, 180, 144, 97} GeV. All UED couplings are 
assumed to be chiral according to Eq. (p^). The leading order cross-section is 0.36 pb 
for SUSY and 2.3 pb for UED. The actual event rate for both sbottoms/KK-bottoms to 
decay according to Eq. ( |4.2| ) highly depends on the decay branching ratios. For UED, 
the KK-bottom has a branching ratio ~ 1/3 to KK-Z^ and KK-Z^ has a branching 
ratio ~ 2/3 to KK-e or KK-jj,, therefore the effective cross-section is 0.11 pb. For 
SPSla, X2 dominantly decays to the stau, and the branching ratio of X2 ~^ yu/e is only 
12%, which makes the effective cross-section much smaller: ~ 5.8 fb. Of course, this 
suppression is not generic, and if necessary one can also consider the stau, though with 
less precision For a similar spectrum (with the 515 GeV sbottom replaced by 

a 565 GeV squark), it was shown in Ref. |26| that the masses can be determined with a 
few GeV uncertainties using 400 events, independent of the spins of the particles. For 
spin determination, the needed number of events is larger. Therefore, we simply use 
the errors given in Ref. in our fit: 



6m N = 4 GeV, 6mx = 4 GeV, Sniy = 4 GeV, Smz = 6 GeV. (4.3) 

We smear the visible particles' momenta according to Appendix 0. About 74% 
(72%) events passed the pr and r] cuts for SUSY (UED). We then apply the likeli- 
hood method described in Appendix ^ to reconstruct the momenta of the two missing 
particles. It allows us to obtain a minimum (Xmin) each combination of the 
visible particles. For simplicity, we only keep the combination that gives the smallest 
Xmin5 although sometimes more than one combinations yield good fits. A fit quality cut 
x"^ < 10 is applied on the events, which further reduces the number of events to 61% 
(60%) for SUSY (UED) with respect to the original number without cuts. Out of the 
final events after all cuts, about 46% events (for both SUSY and UED) have the correct 
combination as checked with the event records from the Monte Carlo simulation. 

After obtaining the missing particles' momenta, we can calculate the momenta of 
the sbottoms/KK-bottoms and examine the production angle, i.e., the angle between 
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Figure 10: The production angle in the center of mass frame of sbottoms/KK-bottoms. 
Left: Monte Carlo truth. Right: from reconstruction. 



and the beam in the center of mass frame. The production angle distributions 
from both the Monte Carlo truth and the reconstruction are shown in Fig. |iy. We 
see that although this distribution is useful to tell the two models apart, the shapes 
of the curves are not as distinct as those of x~^X2 ^^'^ W^'^^'^Z^^^ productions studied 
in Section ^. This is expected. The x~^X2 production processes are 

dominated by s-channel processes from ud initial state. Therefore, it has a simple 
partial wave structure, and the spin of the final state particles in the 2-^2 process 
determines the angular distribution. On the other hand, bb* and b^^^b^^^ productions 
receive contributions from both gg and qq initial states. In addition, gg initiated 
production processes receive contributions from s, t and u channels. Therefore, the 
dependence on the spin of the final state particles in the 2 — )■ 2 process is weakened. 

Additional information can be obtained by studying spin correlation between the 
two decay chains. Of course, there is no correlation in the sbottom pair because they 
are scalar fields. On the other hand, we do expect correlations between the two KK- 
bottoms' helicities. This is analogous to the tt spin correlation ETl except that the 



mass is different. Therefore, we refer readers to Ref. |^ for the detailed discussion. 
The correlation can be observed by examining the opening angle, ^open 0], between 
the two 6-jets. 6'open is defined as follows. We denote unit vector pb as the direction 
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Figure 11: The opening angle between the two b-jets, in the rest frames of the two KK-6's 
respectively. Left: Monte Carlo truth; right: reconstructed. 



of the bottom quark in the rest frame of fe*-^-*, and similarly pi as the direction of h in 
the rest frame of B*-^-*. Obviously ph and pi carry information of the polarizations of iP^^ 
and h'^^\ respectively. One possible observable which can characterize their correlation 
is the product cos6'opcn = Pb ■ Pi- The distribution of the variable can be written as 

= 1 + L'cOS^open, (4.4) 



d cos 6'open 



where D is a constant. The opening angle distribution is shown in Fig. |TT] for both 
Monte Carlo truth and the reconstructed one. As expected, we obtain a flat distribution 
for SUSY and a slope for UED from the Monte Carlo truth, but the reconstructed dis- 
tributions are again modified by experimental cuts and wrong combinations/solutions. 
The constant is D is positive, indicating that the 6-jets tend to move in the same di- 
rection. This can be understood as follows: the KK-bottom pair is more copiously 
produced from the gluon-gluon initial state than from qq. Near the threshold, the final 
state has no orbital angular momentum, therefore has a total angular momentum of 
or 1. The initial gluons do not have orbital angular momentum either and cannot 
form a spin-1 state. Therefore the final state must have total spin 0, corresponding 
to and of the same helicity. Due to the chiral coupling between h and 
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the resulting b and b tend to go in the same direction. That being said, we note that 
the qq initial states contribute significantly to the total cross-section (37%, comparing 
with 15% for tt) and give an opposite slope, which dilutes the effect and makes it a 
difficult measurement. One can also look at the opening angle between the leptons. 
However, unlike the tt case where the charged lepton has the best distinguishing power, 
in our case, the charged lepton gives no advantage due to the fact that the charge of 
the sbottom is not correlated with the charge of the near lepton. 

One can combine the measurements of the production angle and jet-jet opening 
angle to optimize the distinguishing power. For example, we can define a "central- 
forward" asymmetry for the production angle 

(t(| cos6'prod| > 0.5) - o-(| cos6'prod| < 0.5) 

^prod — , 14. OJ 

Ctotal 

and a "forward-backward" asymmetry for the jet-jet opening angle. 

_ (t(cos 6*05011 > 0) - a(cos 6*05011 < 0) 

Ajj — . (^4.bj 

(J total 

The expectation values for the asymmetries after event reconstruction are given by 

^prod(SUSY) = 0.042 (0.078), Aprod(UED) = 0.119 (0.164); 
Aj-j(SUSY) = 0.003 (0.000), Ajj(UED) = 0.026 (0.023), (4.7) 

where the numbers in the parentheses are from the Monte Carlo truth for comparison. 
The asymmetries are small, in which case the statistic errors are simply given by 1/y/N 
where is the number of available events after cuts. By combining the two measure- 
ments, it is possible to distinguish the two spins with ~ 1700 events at 95% level after 
cuts^. This number is obtained by assuming a sample of pure signal events. In reality, 
more events may be required due to systematic uncertainties, the SM backgrounds, 
as well as contaminations from other new physics processes with the same final state 
particles. Of course, one should combine other information from the decay chain to 
better determine the models. For example, the jet-near lepton invariant mass gives us 
information about the spin of the Xg/KK-Z^. For the distribution is fiat since we 
do not know the charge of the 6-jets, while for KK-Z^ a second order polynomial can 



be seen P, |I2|. However, we emphasize again that this kind of information is not a 



direct measurement of the spin of the first particle in the decay chain. 

^Ignoring a small correlation between Aprod and Ajj, we define x^(SUSY) = 
[Ap,od(exp) - Ap,od(SUSY)]Va2^„d + [Ajj(exp) - Ajj(SUSY)]V^Tfj and x'(UED) = [Aprod(exp) - 
Aprod(UED)]^/crpj.Qjj + [Ajj(exp) — Ajj(UED)]^/CTj?j, where Aprod(exp) and Ajj(exp) are the experimental 
values and dprod — crjj = 1/\/N. We estimate the needed number of events by requiring that the 
correct theory has 95% probability of having the smaller and hence being selected. 
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5. Discussion and Conclusions 

In many scenarios of new physics beyond the Standard Model, such as supersymmetry 
and UED, the decay of new physics particles frequently leads to long decay chains 
ending with a stable massive neutral particle with undetectable momentum. In the 
article, we have studied methods of reconstructing the kinematics of such decay chains. 
We began with the assumption that the masses of the new physics particles involved 
in the decay chain have been measured. We then showed that their momenta, in 
particular the momentum of the stable neutral particle, can be fully reconstructed. As 
an application of this method, we used the kinematic information to determine the spin 
of new particles and showed that different new physics scenarios, such as supersymmetry 
and UED, can be distinguished with this method. Well studied methods using Lorentz 
invariant variables are not directly applicable to the measurement of the spin of the 
particle at the first or the last step of the decay chain. With full kinematic information 
of the decay chain, we are able to probe the spin of the particle which initiates the 
decay chain. 

We performed two case studies. First, wc considered the kinematic reconstruction 
and spin measurement with information from only one side of the event, i.e., one 
decay chain. We also demonstrated a "double chain" analysis, using the kinematic 
information to obtain the production angle and the correlation between the decay 
products of two new physics particles, one on each side of the events. We expect 
these two methods to be complementary. The final state in the single chain analysis is 
obvious more inclusive. At the same time, for the distribution of the decay products to 
contain useful spin information, the particle under consideration needs to be produced 
in a polarized state, and its coupling to its decay product has to be chiral. On the 
other hand, extraction of non-trivial spin correlation in the double chain analysis can 
be successful without such special requirements on the couplings of new physics states. 
However, it obviously requires precise knowledge of both sides of the decay chains. 

We have demonstrated our method using a set of particular benchmark models, 
SPSla with a light gravitino and a model with similar mass spectrum in the cases 
of supersymmetry and UED, respectively. While these benchmark models are not 
designed to allow an easy spin measurement, our choices of the production channel and 
the strategy do take advantage of specific features of the spectrum. For example, we 
have relied on the fact that squark/KK-quark exchange does not contribute significantly 
to the X2X^/^*^^'*^^^^^ production. Similar method should be applicable to other 
models, even though the specific choice of channels and strategy can be different. We 
remark that this situation is expected to be quite generic. Due to subtleties in the 
extraction of spin information and the virtually infinite number of possibilities of new 



- 20 - 



physics models, it is impossible to have an observable which is universally applicable. 
However, once enough details of the new physics states, such as masses and quantum 
numbers, are known, it is likely that specific variations of several proposed classes of 
spin measurement methods, such as the method demonstrated here and the invariant 
mass method, can be adapted to accomplish the task. 

To focus on our demonstration of the reconstruction and spin measurement method, 
we have used exclusive samples of signal events in our analysis. Realistically, achiev- 
ing an exclusive sample of high purity requires strict cuts to suppress the Standard 
Model backgrounds, and contaminations from other new physics channels. Performing 
a careful study of the reach in the specific examples discussed here could be interesting. 
However, the design and optimization of such cuts will inevitably be very model depen- 
dent. As a result, the conclusion of such a study is less likely to be representative in a 
large class of models. Therefore, such studies will be more effective after particular new 
physics channels have been identified at the LHC. Due to the expected low efficiency in 
isolating such exclusive samples, and the prerequisite of mass measurements, we expect 
the method presented here will be useful only with large statistics. 

Instead of assuming prior knowledge of the masses, we could in principle perform 
a combined strategy which fits both masses and spins. Moreover, reconstruction of 
momenta can also help measure other properties of the new physics, such as the chirality 
of the couplings as we have already alluded to in Section ^ We will postpone further 
development of our methods in those directions to future studies. Finally, we note 
that momentum reconstruction is also useful for observing CP violations It is 

interesting to study possible applications of the methods presented in this article. 
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A. Minimization Method for Momentum Reconstruction of 
Over-constrained Systems 

We describe the general formulation of momentum reconstruction for over-constrained 
systems such as the three-step cascade decay chains depicted in Fig. |[ We also point 
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out several subtleties that arise from our assumptions about the measurement variables 
and address how to resolve these subtleties. 

The goal of event reconstruction is to determine the invisible particle momenta 
that maximize the likelihood (by minimizing x^) of hypothetical mass shell relations of 
the involved particles. The ingredients can be summarized as follows: 

• Parameters: The parameters are the quantities to be reconstructed. They are 
varied in the reconstruction procedure. Here, the parameters are the momenta of 
the invisible particles, which we denote as 9i, I — 1,. . . ,N. 

• Measurement variables: The measurement variables are fixed quantities for each 
event, such as the momenta of the visible particles (jets and leptons). We denote 
these quantities by x^. A; = 1, . . . , m. 

• Nuisance parameters: These quantities are shared by all events and are assumed 
to be known a priori. Here, the nuisance parameters are the mass parameters 
of the particles involved in the event topology. For simplicity, in our analysis we 
will treat these parameters as measurement variables {i.e., as part of the x^'s). 

• Hypotheses: The hypotheses are the relations that should be satisfied by the 
parameters, measurement variables, and nuisance parameters of the system. Here, 
these are the mass-shell relations and the missing momentum sum. We will label 
the hypotheses as yi(xfc; 9i) = 0, i = 1, . . . ,n. 

The procedure is to calculate the likelihood of the hypotheses and maximize the 
likelihood with respect to the parameters. Assuming Gaussian statistics, the maximum 
likelihood condition can be equivalently obtained by minimizing of th^ hypotheses 
for a given event. Under this assumption, the statistical dependence between the 
measurement variables can be described by the covariance matrix: 



Here (q) and q denote the statistical mean value of some arbitrary variable q, evalu- 
ated over statistical ensemble of experimental measurements with the same physical 
configuration. 

To define the of the hypotheses {yj = 0}, one defines the covariance matrix Vij 
of Yi at a given (x^, 6i) by 



Ukl = ((Xfc - Xfc)(x; - X;)). 



(A.l) 



Vij = {yiYj) 



(A.2) 




(A.3) 
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where we have used the fact that y^'s must have zero mean value. Since the covariance 
matrix Vij is not diagonal in the systems under consideration, the function of our 
hypotheses is given by 

x'{9j) = [yi^,9)r.V-'-[yi^,9)], (A.4) 

where [y] denotes a column vector constructed from the yj. 

For the cascade decays depicted in Fig. |, the kinematic constraints of the system 
can be summarized by the following set of equations: 



yi 


= pj - m% = 0, 




(A.5) 


y2 


= (Pi +P3f -rux = 0, 




(A.6) 


ys 


= {Pi +P3+P5)'^ - my 


= 0, 


(A.7) 


y4 


= (Pl +P3+P5+P7f - 


m| = 0, 


(A.8) 


y5 


= pI - m% = 0, 




(A.9) 


ye 


= {.P2+Pif - = 0, 




(A.IO) 


y? 


= {P2 +P4:+P&f - my 


= 0, 


(A.ll) 


ys 


= {P2 + Pi + P6 + Psf - 


m| = 0, 


(A.12) 


yg 


= Pt+P2-f =0, 




(A. 13) 


yio 


= p\ + pI — = 0. 




(A. 14) 



For this system, the parameters are the momentum variables for the invisible particles 
in each event, which are the eight real quantities p^ and P2 {fi = 0, . . . ,3). There are 30 
measurement variables for each event, which include the visible particle momenta 
Ps y P5 y P7 y P4t Pe T Ps y missing transverse momentum '^^ and and the 4 "nuisance 
parameters": rriN, mx, my and rriz- The visible particle momenta are measured in 
terms of the transverse energy Et = A/pl"+~p|~+m?, the pseudo-rapidity t] and the 
azimuthal angle 0. Since the visible particles in this case are massless, we can write 

pf" = {ErCr,, Etc^, Ets^, Ets^), (A.15) 

where c<^ = coscp, = sincp, Cr, = coshr/, and Sn = sinhr/. Denoting the statistical 
errors by ^^j,, 5^ and 5^ for Et, r] and (p, respectively, the covariance matrix is given by 

{P,P.) = (A. 16) 



-23- 



Electrons and muons [53, H : 


SpT_ _ 
PT 

69 = 


< 2.4, Pt > 10, 
= 0.008 ©O.OOOlSpT, 
= 0.001,50 = 0.001. 


Photons [ 


44, 


45 




\r] 

5E 
E 

5rj = 


< 3.0, Pt > 10, 
^©^©0.0026, 
0.0011, 50 = 0.003. 


Jets: 


61] = 


\ri\ < 3.0, PT > 100 GeV, 
f M © ^ © 0.033, for \r]\ < 1.4, 
1^ ^8 © ^ © 0.043, for 1.4 < \r]\ < 3.0, 
7 = 0.03, 50 = 0.02 for \ri\ < 1.4, 
0.02, 50 = 0.01 for 1.4 < \ri\ < 3.0. 



Table 2: Experimental errors from the measurements of particle momenta: In our analysis, 
parton level events are smeared according to the above Gaussian errors. The observables 
of energy dimension are in GeV units and the angular and the rapidity variables are in 
radians. Simple acceptance cuts on pseudo-rapidity r] are also applied. For electrons and 
muons, the resolution here corresponds roughly to the CMS tracking system performance in 
the central region (r/ = 0) |^3|. The resolution becomes slightly worse at higher rapidity until 
I'?! ^ 2 where it starts to diverge. We ignore this rapidity dependent effect. For photons, the 
resolution for position measurement corresponds to the CMS ECAL performance obtained 



using electron beams (10 < < 50 GeV) in Ref. [45| 



Here we have omitted the particle index for the same particle. For different particles 
i and j we have (p^p^) = 0, because the momentum measurements are uncorrelated 
for different particles. On the other hand, the missing transverse momentum measure- 
ment is correlated with the other measurements. The entries of the covariance matrix 
involving and are thus given by 



(A.17) 



where fi, v in Eqs ( A.17| ) are restricted to x and y. The experimental errors from the 
measurements of particle momenta and missing transverse momentum are summarized 
in Table 

With these ingredients, the next step is to determine as given in Eq. ( |A.4[) . 
However, one immediately encounters a difficulty in doing so because the covariant 
matrix V is singular and the from Eq. ( |A.4D is ill-defined. The physical reason for 
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this singularity is that the two decay chains in Fig. ^ are assumed to be symmetric 
and the masses of the particles in both chains are identical. In a realistic situation, 
the assumption of symmetric chains can be incorrect due to effects such as finite decay 
widths and particle misidentification. Nevertheless, it is desirable to have a method to 
resolve this problem without abandoning the symmetric chain assumption. 

Our procedure for addressing the situation is to introduce a regulator that controls 
the divergence as follows. We double the number of nuisance parameters by assuming 
that the masses of the particles N, X, Y, Z in each chain are effectively independent 
variables, which leads to a covariance matrix e.g. for rriN of the form 

^2' 11 



{{TTiNi - mNi){mNj - niNj)) = (Sttin) j ■ (A.18) 

We then apply a very small perturbation e to Eq. ( |A.18| ), 

{{mm - rnm){mNj - rn^)) = {SttlnY' {^^^_^ ■ (A.19) 

After this perturbation, Vij is no longer singular and is well-defined. The numerical 
results presented in Section || are obtained by setting e = 0.01. The divergent part of 
can be extracted by taking e — )■ 0, 

x\e) = xm + -m- (A.20) 

The regular part of x^(^), which has been denoted as x^{(^)i can be extracted from the 
pseudoinverse of the symmetric matrix Vij. The pseudoinverse V is defined by 

V = W-D-W-^, (A.21) 

where W is the diagonalizing matrix of V: 

V = W -D- W~\ with D = diag(rfi, c/2 • • • c^m, 0, 0, . . . 0), {di ^ 0), (A.22) 

and D is given by 

D = diag{d^\d^\...,d;;^\0,..). (A.23) 
The regular part of is then given by 

7' = [yf-V-[y]. (A.24) 
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The 1/e term plays the role of a penalty term since it effectively confines the configu- 
ration space (parameter space) to the solution space of h{9) = 0. To obtain this term, 
we define the projection matrix P of Vij onto the space of zero eigenvalues as follows: 

P{V) = W ■ diag(0, 0, ....0, 1, 1, 1, ..) ■ W-\ (A.25) 

where the I's correspond to the (n — m)-dimensional subspace of V with zero eigenval- 
ues. h{9) in Eq. ( |A.2CI|) is then given by 

n / n \ ^ 

h{9) = [yf-PiV)-[y]= [J2ywJ . (A.26) 

k=m+l \i=l / 

In summary, the minimization must be carried out in the restricted space of 

n 

5^y«W^^fc = 0. (A.27) 

i=l 
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